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Abstract 

The author's recent results on an asymptotic description of the Schlesinger equa- 
tion are generalized to the JMMS equation. As in the case of the Schlesinger equa- 
tion, the JMMS equation is reformulated to include a small parameter e. By the 
method of multiscale analysis, the isomonodromic problem is approximated by slow 
modulations of an isospectral problem. A modulation equation of this slow dynam- 
ics is proposed, and shown to possess a number of properties similar to the Seiberg- 
Witten solutions of low energy supersymmetric gauge theories. 



KUCP-0106 



*Partly supported by the Grant-in- Aid for Scientific Researches, Priority Area 23f "Infinite Analysis" , 
the Ministry of Education, Science and Culture, Japan 



1 



1 Introduction 



This paper aims to supplement the author's recent work on an asymptotic description of 
isomonodromic problems |l[ in the language of Whitham equations p|. In the previous 
paper, the author considered the Schlesinger equation ||, and proposed a modulation 
equation that is expected to govern the Whitham dynamics of an underlying isospectral 
problem. The modulation equation and its solutions exhibit remarkable similarity with the 
so called "Seiberg- Witten solutions" of low energy supersymmetric gauge theories |4] . In 
this paper, we shall extend these results to the Jimbo-Miwa-Mori-Sato (JMMS) equation 
||. The JMMS equation is the simplest generalization of the Schlesinger equation with 
an extra irregular singular point at infinity, and has found many applications in solvable 
lattice models, field theories, and random matrices |5|, || [7]. Furthermore, this equation 
possesses a kind of duality that exchanges the roles of regular singularities at finite points 
and the irregularities at infinity |§. This duality is an isomonodromic counterpart of a 
similar duality in isospectral problems of Gamier and Moser type || . A main concern of 
the following consideration is the fate of this duality in the modulation equation. 



2 JMMS Equation 

In this section, we collect basic facts about the JMMS equation from the literature f|, |7], || . 
The JMMS equation governs isomonodromic deformations of the matrix ODE system 

^- = M(X)Y (1) 

with an r x r matrix coefficient 

N * m 

M(A) = tf + ETzV' (2) 

i=l % 

where A{S are r x r complex matrices, and Ui a diagonal matrix of the form 

r 

U = diag(«i, • • -,u r ) = X] u a E a . (3) 

a=l 

(Here E a is the matrix whose 7) element is 8 a p8 ai .) U and u a are "time variables" 
of isomonodromic deformations. A = tj (i = 1, • ■ ■ , N) are regular singular points, and 
A = 00 an irregular singular point of Poincare rank 1. In order to avoid logarithmic terms, 
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we assume, throughout this paper, that the eigenvalues 6i a (a — 1, • • • , r) of A* have no 
integer difference, and that u^s are pairwise distinct: 



Under this assumption, isomonodromic deformations are generated by the PDE's 



(4) 



dY 



' Y, = (XE a + B a )Y, 



X-t 



where 



N 



(5) 



(6) 



0&a) Ua ~ U P i=l 

The Frobenius integrability conditions give the so called "zero-curvature" (or "Lax") 
representation of isomonodromic deformations. The essential part of the zero- curvature 
representation consists of the equations 



d A; d , 

— h —, 7TT- — M(X) 



o, 



— XE a — B a , — - M(A) 
ou a oX 



dU X — U dX 
They give the following equations for A^s (JMMS equation): 



= 0. 



(7) 



dAj 
~dt~ 

dA^ 



[tjE a + B a , Aj 



ti tl 



The isomonodromic deformations have two sets of invariants. One of them consists of 
the eigenvalues 9 ja of A/s: 



= 0. 



(9) 



(10) 



dti du c „ 

Another set of invariants are the diagonal elements of A^: 

dU du Q 

These invariants are nothing but characteristic exponents of the ODE at the singular 
points. 

The invariance of 6 ia is related to a coadjoint orbit structure behind the JMMS equa- 
tion. Note that the right hand side is a commutator of the form [• • • , Aj]. This implies 
that the motion of Aj is constrainted to a coadjoint orbit Oj in gl(r, C); the eigenvalues 



6j a are invariants that determine this coajoint orbit. Thus, the JMMS equation is a non- 
autonomous dynamical system on the direct product 0\ X • • • x of coajoint orbits. 
This dynamical system can be written in the Hamiltonian form 

f)A ■ 8 A ■ 

with the Hamiltonians 

^ = Tr (uAi + £ AlA] 



j(^i) ti tj , 

K a = Tr( E a jt,UA i + £ EaA °° E P A A . (12) 
The Poisson bracket is the standard Kostant-Kirillov bracket on gl(r, C): 



3 Dual Formalism 

In this section, we review the notion of dual isomonodromic problem ||. Although we 
have assumed the genericity condition on characteristic exponents, the dual problem itself 
can be formulated in a more general setting. (This is also the case for isospectral problems 
0.) Suppose now that the rank of A4 takes a general value £j. is also a coajoint orbit 
invariant, and constant under the JMMS equation in the generalized sense. 

The first step towards the dual problem is to express the rational matrix M(A) as: 

M(A) — U- t G(\I — T) -1 F, (14) 

where F = (F aa ) and G = (G aa ) are t x r matrices, i = an d T is a diagonal matrix 
of the form 

N 

T = E^A, A = H h (15) 

8=1 

In particular, Aj is written Ai = — l GDiF. F aa and G aa may be understood as canonical 
coordinates with the Poisson bracket 

{F aa , Fb/s} = \G aa , G\,p} = 0, {-F aQ , G^} = 5 a b5 a p. (16) 
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The map (F, G) \— > (Ax, ■ • • ,An) then becomes a Poisson map, and the JMMS equation 
can be derived from a Hamiltonian system in the (F, G)-space with the same Hamiltonians 
Hi and K a . 

The dual problem is formulated in terms of the rational i x I matrix 

L(n) = T - G(fj.I - U)- U G, (17) 

which can be rewritten 

m=T+j:-^- as) 

a H 1 "'a 

with P a = —FEjG. Note that P a is a rank 1 matrix; this restriction, too, can be relaxed 
0, |9|. The map (F, G) i— > (Pi, P2, • ■ ■) is a Poisson map for the Kostant-Kirillov bracket 

{P a ,ab, Pf3,cd} = 5 a f3( — 5bcP/3,ad + Sd a P/3,cb)- (19) 

The previous Hamiltonians Hi and K a can also be viewed as a function of Pa's, and 
induces a non-autonomous Hamiltonian system in these new variables. This again gives 
an equation of JMMS type. The corresponding isomonodromic problem can be written 

% = -Lfr)Z, ¥-=-*^Z, ^ = -(^D i + Q i )Z, (20) 
a A ou a /i — u a oti 

where 

Qi = ~E DiP °° D i+fi P °° D \ P oa = -'£P a . (21) 

4 Introducing Small Parameter 

We now introduce a small parameter e by the following substitution rule (which was first 
proposed by Vereschagin for the Painleve equations |10| ) : 



d d d d d d 



d\ <9A' dti 9Tj' du a du 
1 1 



x _ f X _ T - E^^E^- ( 22 ) 

As we shall show below, T« and U a play the role of "slow variables" in the terminology of 
multiscale analysis. The JMMS equation now takes the e-dependent form 



dAi [AijAj^ 



- = (I —O,;) • = •(>,, U + ^ ,Aa 

dA 



[TjE a + B a ,Aj] , (23) 



dU a 
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where 



(24) 



U a -U p 

A central idea of multiscale analysis is the use of two sets of independent variables 
("fast" and "slow" variables) that represent two different scales of phenomena. In the 
present setting, A/s are thus first assumed to be a function of "fast variables" (t,u) = 
(ti,u a ) and "slow variables" (T, U) = (T i: U a ), then restricted to the subspace 



(25) 



Accordingly, derivatives in the original differential equation are replaced by the sum of 
derivatives in the fast and slow variables: 



dA, OA, 



dTi 



+ e 



dA, 



OA, 



dT,. ' dU a 



du„ 



+ e 



dAj 



(26) 



Subsequent analysis is a kind of perturbation theory Namely, one assumes the asymp- 
totic form (as e — > 0) 



Aj = Af (t, u, T, U) + Af (t, u, T,U)e + 



(27) 



and considers the differential equation order-by-order in e-expansion. The lowest order of 
this expansion for the JMMS equation yields the following equation: 



8A 



(0) 



dti 

dAf 
dU„ 



= a 



l3> Ti-Ti 



|(0) 



fc(#i) 



Ti-Ti 



(0) 



t^ + b^Ka^ 



(28) 



Here is given by the same formula as B a except that Aoo is replaced by A^ = 

Note that the lowest order equation itself does not determine the (T, ?7)-dependence 
(i.e., "slow dynamics") of Af^s. The slow variables appear as parameters in the lowest 
order equation. In the standard recipe of multiscale analysis, the slow dynamics is to be 
determined by the next order analysis. We shall, however, take a different approach to 
bypass technical difficulties (or, rather, complexity). 
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5 Isospectral Problem 



Our approach to slow dynamics, which is still heuristic but considerably transparent, is 
based on the fact that the above lowest order equation is an isospectral problem. The 
above equation indeed has the isospectral Lax representation 



d_ 4 (0) 

dU + A - Ti 



M (0) (A) 



0. 



_d_ 

dUr 



\E a -BW,M<®(\) 



0. 



where 



N 



M (0) (A) = U + J2 



A 



(o) 



(29) 



(30) 



In particular, the characteristic polynomial detf/il— M^(A)J is constant under the (t, u)- 
flows. It however depends on (T, U). Thus the slow dynamics in (T, U) may be described 
as slow deformations of the characteristic polynomial or of the spectral curve 

det(/i/-M (0) (A)) = 0. (31) 

Spectral curves of the above type has been studied in detail in the context of isospectral 
problems of Gamier and Moser Jll], |i~2"f . A fundamental result is that the spectral curve 
Co on the (A, //)-plane can be compactified to a nonsingular curve C. This is due to the 
previous assumption on the eigenvalues of Aj's and U ; if they have multiple eigenvalues, C 
has nodal singularities. Let 7r denote the projection 7r(A, /i) = A from Cq to the punctured 
Riemann sphere CP 1 \ \T\, • • • , Tjv, oo}. This extends to C and gives an r-fold ramified 
covering of C over CP 1 . C is obtained from Cq by adding several points to the holes over 
A = Ti, • • • , T^, oo. The holes over A = oo are filled by the r points (A, /i) = (oo, U a ) 
(a = 1, ■ • ■ , r). The situation of holes over A = Tj becomes clearer in terms of the new 
coordinate 



TV 



fi = f(\)fM, /(A) = n(A-T,), 



i=i 



in place of fi. The equation of the spectral curve can be rewritten 



F(A,/2) =det(/LT-/(A)M (0) (A) 



0. 



(32) 



(33) 



Compactification is achieved by adding the r points (A, fl) = (Tj, f'(Ti)6 ia ) (a = 1, • ■ • , r) 
The compactified spectral curve has the genus 



g=\{r-l){rN-2). 



(34) 



In order to describe slow dynamics, one has to choose a suitable system of parameters 
( "moduli" ) in the space ( "moduli space" ) of all possible spectral curves. Slow dynamics 
can be formulated in the form of differential equations ( "modulation equations" ) for those 
moduli with respect to the slow variables (T, U). Note that the isomonodromic invariants 
9 ia and A^^ are also constants of slow dynamics. 

Such a set of moduli can be singled out from the modified characteristic polynomial 
F(X, p). F(X, p) turns out to have an expansion of the form 

F(X, p) = F°(A, p) + /(A) £ E* h ms \ m p\ (35) 

s=2 m=0 

where 5 S is given by 

8 S = (N - l)s - N. (36) 

F°(X,p) is a polynomial whose coefficients are determined by the isomonodromic invari- 
ants and (T,U) only. The coefficients hs s -i s (2 < s < r), too, are comprised of these 
quantities only. The remaining coefficients h ms (2 < s < r, < m < 5 S — 2) give the 
moduli in the above sense. Thus, as in the case of moduli spaces of Seiberg-Witten type 
||, we have altogether g moduli: 

I>« - 1) = 9- (37) 

s=2 

Thus, the leading order term of the multiscale expansion turns out to be a solution 
of the isospectral problem with slowly varying spectral invariants h ms = h ms (T, U). Since 
this kind of isospectral problems can be solved by Abelian (or Riemann theta) functions 



TT| , the multiscale expansion A± = Af^ + 0(e) gives an asymptotic expression of 



JMMS transcendents as a modulated Abelian functions. 



6 Modulation Equation 

The first step for deriving the modulation equation is the following WKB ansatz (which 
dates back to Flaschka and Newell's early work on isomonodromic problems flT3"|): 

Y = [t, u, T, U, A) + (1) (*, u, T, U, A)e + • • •) exp e^^T, U, A) . (38) 
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Here Y and 0^'s are understood to be vector valued; S is a scalar phase function. Note 
that whereas depend on both the fast and slow variables, S is assumed to be a 

function of slow variables only. Inserting the above ansatz into the isomonodromic linear 
problem and expanding in powers of e give a series of differential equations for 0( fc )'s and 
S. 



The lowest order equation reads: 

dS 
dX 

<90 (o) dS_ 
~dt~ + dT i 
<90 (o) dS 
du a dU a 

(Note that these equations determine 
function h = h(T, U) of slow variables 
isospectral linear problem 



M(°)(A)0W 
dA (o) 



A — Tj 



(39) 



up to multiplication (p^ i— ► by a scalar 
These equations can be rewritten into the 



^ = M (0) (A)^, 
by change of variables as 

dS 

^ = SA ' 



#0 



A- 



(0) 



X — Ti 



■1>, 



<9-0 

dUrv 



OS 



as 



(40) 



(41) 



In particular, the characteristic equation det^u/— M^°\X)) = is satisfied by = dS/dX. 

We derive the modulation equation by identifying this with the algebro-geometric 
expression of ip in terms of Baker-Akhiezer functions (with a suitable correction by a 
scalar factor h such that = (j^°>h). In the present case, such an expression takes the 



form [11 



i) = 0exp(^^ Ti +J2u a Q Ua ), (42) 

where is a vector valued function made from various Riemann theta functions on the 
spectral curve, and Qx t and Qu a are primitive functions of meromorphic differentials dVt^ 
and Q Uoi , 

/-(A,m) 

Qt. = J dn Ti , n Ua = J dn Ua . (43) 

[More precisely, such an expression of Baker-Akhiezer functions is available upon fixing 
a symplectic basis Ai, Bj (I = 1, ■ ■ ■ , g) on the spectral curve. In the present situation, 



the symplectic basis has to be chosen so that the above WKB solution gives a correct 
approximation to the isospectral problem; this is a nontrivial problem. In the following, 
we assume that such a suitable symplectic basis is selected.] The meromorphic differentials 
dVtTi and dQu a are characterized by the following conditions: 

• dflTi has poles at the r points in 7r~ 1 (T i ) and is homomorphic outside 7r _1 (Tj) . Its 
singular behavior at these points is such that 



dQ^ = —dfi + non-singular. 



(44) 



dVLu a has poles at the r points in tt 1 (oo) and is holomorphic outside tc 1 (oo) . Its 
singular behavior at these points is such that 



dQ Ua = d\ + non-singular. 



(45) 



dQ^ and dQjj a have vanishing Aj-periods: 



dfli 



dVLj 



(/ 



g)- 



(46) 



Matching the exponential part in the two expressions of xp, we obtain the equations 



dS_ dS_ 
W l ~ T " dU~ a 

This is an expression of our modulation equation. 

This modulation equation can be rewritten 



VLrr. 



d_ 

df, 



dS 



dQ 



d 



Ti 



A=const 

in terms of the meromorphic differential 



dU a 



dS 



dQ 



A=const 



dS = fld\. 



(47) 



(4£ 



(49) 



Here, as in the literature of Whitham equations fl~4] , pT5| , "A = const " stands for differen- 
tiation while leaving A constant. In particular, the following equation of Flaschka- Forest- 



McLaughlin type [16] are satisfied under the above equations: 



_d_ 

df, 



dQn 



A=const 



S—dVtjj 

dTi Ua 



A=const 



A=const 







dU r 



dVL 



Ti 



A=const 



A=const 



_d_ 



dtt 



A=const 



(50) 
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Thus our modulation equation turns out to give a special family of Whitham equations. 
These equations, though rather indirectly, determine the slow dynamics of the spectral 
curve (hence of the spectral invariants h ms = h ms (T, U)) for the isospectral aproximation 
Ai ~ Af ] to the JMMS equation. 

The same multiscale argument applies to the dual isomonodromic problem, too. The 
WKB ansatz now takes the form 

Z = ( x (°) (t, it, T, U, //) + (t, u, T, U, y)e + ■ ■ •) exp e~ l E(T, U, /i). (51) 

This leads to slow dynamics of the dual isospectral problem |§ based on the dual expres- 
sion 

det ( XI - L^(/4) = (52) 



dT,, 



= dn Ua , (53) 

/i=const 



of the same spectral curve. The modulation equation is obtained in the dual form 

= dVtr , d 
dU a 

where 

rfS = -Xdfx. (54) 

It should be noted that a similar dual representation is known for the Seiberg-Witten 
solutions 0. 

7 Period Map and Prepotential 

The above modulation equation, like the modulation equation of the Schlesinger equation 
|J, can be treated as a variant of problems of Seiberg-Witten type. We show a list of 
main results below. 

1. General solutions of the modulation equation are obtained by an inverse period 
map. Let us consider the ^-dimensional period map h = (h ms ) i— > a = (aj) given by the 
^/-periods 

ai= I dS (I = !,•••,<?). (55) 



J A x 

On very general grounds |l] , one can prove that the Jacobian of this map is non- vanishing. 
The slow variables (T, U) are included in the period integrals as parameters. Accordingly, 
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the the inverse period map should be written h = h(T,U,a). With a being fixed, this 
gives a family of deformations of the spectral curve, which turns out to solve the above 
modulation equation. Varying a, one obtains a general solution. 

2. With (T, U) being fixed, this inverse period map yields deformations of Seiberg- 
Witten type. This means that the meromorphic differential dS obeys the deformation 
equation 

d dS =doo I (/ = !,- ••,<?), (56) 



dai 

A=const 

where dujis are a basis of holomorphic differentials on the spectral curve normalized as 



f dcoi = 5ij. (57) 

JAj 



3. A prepotential of Seiberg-Witten type exists. The prepotential T = J r (T,U,a) can 
be defined as a (non-zero) solution of the following equations: 

da~r hu df= H - m: =Ka - (58) 

Here 6/ denotes the ^/-period of dS, 



h = j B dS, (59) 



and Hi and K a are the Hamiltonians already introduced. In particular, the second deriva- 
tives of T coincide with the matrix elements of the period matrix: 

d 2 T 



da i da j 



= r u = 



I dui. (60) 
Jb.t 



8 Conclusion 



We have shown that the previous work on the Schlesinger equation can be generalized to 
the JMMS equation. Although the derivation of the modulation equation is still heuristic, 
the modulation equation has turned out to possess a number of remarkable properties that 
strongly suggest that our modulation equation is a correct one. We have also seen that 
the duality structure of the JMMS equation is inherited to the modulation equation. 

A straightforward generalization of the present situation is to add one more irregular 
singular point at, say, A = 0. The rational matrix M(A) then takes the form 

M(A) = (/ + E I ^ + ^ 2 . (61) 
12 



In this kind of isomonodromic problems, U and V are usually required to be semi-simple. 
The beautiful duality in the JMMS equation, however, then ceases to exist. Relaxing this 
requirement is an interesting issue. A similar situation takes place in the iV-periodic Toda 
chain; it has an N x N Lax matrix of the form L = U + W\~ l + VX~ 2 with U and V 
being nilpotent. It is well known that a dual representation of this case is provided by 
a 2 x 2 monodromy matrix M on the lattice. This fact seems to suggest that a similar 
duality can exist in isomonodromic problems with two irregular singular points, too. 

Another intriguing issue is to extend the present consideration to the case of an elliptic 
matrix. Remarkably, our M-matrix has almost the same structure as the Lax operator 



of the XXX Gaudin model [17|. Since an elliptic analogue of the XXX Gaudin model is 



given by the XYZ Gaudin model |y| , It is natural to expect that there will be an elliptic 
version of the Schlesinger or JMMS equation. In the case of 2 x 2 matrix, such an elliptic 
isomonodromic problem will be equivalent to Okamoto's isomonodromic problem on a 



torus |TjJ. The corresponding isospectral problem and its Whitham modulations should 
be described by a spectral cover of the torus. This situation, too, is interesting in the 
context of the Seiberg-Witten solutions, because the Seiberg-Witten solution of N = 4 
super symmetric SU(2) gauge theory is also related to a spectral cover of the torus ||. 
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